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Abstract 

PEM-fuel  cells  operated  with  hydrogen  and  air  offer  promising  possibilities  for  the  decentralised  energy  supply  in  stationary  and  mobile 
applications.  But,  there  is  still  a  remarkable  need  for  more  research  for  the  optimisation  of  the  single  components  even  though  the  research, 
especially  the  development  of  membrane-electrode  units  and  bipolar  plates,  has  made  considerable  progress  within  recent  years.  This  also 
applies  to  the  definition  of  suitable  test  algorithms  and  parameters  for  recording  the  characteristics  as  well  as  long  time  tests. 

The  investigations  for  single  cells  or  stacks  can  thereby  be  subdivided  into  investigations  in  the  stationary  and  the  dynamic  state.  This 
paper  shows  a  simplified  approach  for  the  evaluation  and  modelling  of  PEM-fuel  cell  stacks  for  the  stationary  state.  Based  on  the  definition  of 
regression  approaches  for  the  dependence  of  several  parameters  the  number  of  stack  parameters  is  initially  reduced  to  a  sufficient  number  of 
values.  The  remaining  parameters  are  used  to  form  an  energy  model  which  can  be  combined  with  the  energy  models  of  auxiliary  components 
like  air  compressors  and  sinus  inverters  [1]. 

Algorithms  for  recording  the  parameters  still  have  to  be  defined.  Therefore,  the  tests  for  the  preparation  of  a  set  of  characteristic  curves  can 
be  very  time  consuming  because  of  the  multitude  of  variable  parameters.  This  paper  presents  and  discusses  optimised  measurement  algorithms 
for  the  evaluation  of  PEM-fuel  cell  stacks  to  reduce  this  time  exposure. 

©  2005  Elsevier  B.V.  All  rights  reserved. 

Keywords :  Fuel  cell;  Measurement  parameter;  Measurement  algorithm;  Energy  model 


1.  Introduction 

A  fuel  cell  facility  for  stationary  power  generation  consists 
of  a  multitude  of  components  in  addition  to  the  stack.  These 
components  are:  reformers,  if  fuel  processing  is  necessary,  as 
well  as  control  sinus  inverters,  pumps  and  compressors  for 
transporting  the  process  media.  An  optimal  adaptation  of  the 
operating  points  of  the  several  components  to  the  fuel  cell  is 
extremely  important  for  fuel  cell  facilities  in  the  lower  power 
range. 

The  operational  behaviour  of  the  several  components  is 
important  for  adaptation.  For  evaluation  of  the  operational 
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behaviour  of  PEM-fuel  cells  numerous  models  exist  which 
more  or  less  satisfactorily  describe  the  behaviour  of  the  fuel 
cell.  These  models  normally  require  a  lot  of  characteristic 
variables  which  are  not  known  if  a  commercial  stack  is 
used.  For  first  considerations  energy  models  are  sufficient 
which  do  not  require  a  detailed  knowledge  about  the 
investigated  fuel  cell  stack.  However,  the  energy  balance 
of  a  fuel  cell  is  characterised  by  a  larger  number  of  vari¬ 
ables  as  opposed  to  peripheral  components  like  pumps 
or  compressors  which  can  be  described  with  only  a  few 
parameters. 

To  describe  the  fuel  cell  energy  balance  the  parameters 
of  the  process  media  on  the  inlet  and  the  outlet  of  the  stack 
as  well  as  the  voltage-current  behaviour  dependent  on  the 
parameters  of  the  process  media  has  to  be  known.  The  set  of 
parameters  can  still  be  simplified  by  making  some  assump¬ 
tions  without  remarkable  information  losses  regarding  the 
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Nomenclature 

A 

area  (cm2) 

b 

Tafel  parameter  (mV) 

cp 

specific  heat  capacity  at  constant  pressure 
(Jkg-'K-1) 

Cs 

radiation  of  the  black  radiator 

(5.77  Wm“2  K"4) 

Cw 

specific  heat  capacity  of  water  in  the  liquid 
state  (Jkg-1  K-1) 

d 

characteristic  diameter 

F 

Faraday  constant  (96,485  C  mol-1) 

g 

acceleration  of  gravity  (9.81  m  s-2) 

Gr 

Grashof  number  (-) 

h 

specific  enthalpy  (Jkg-1) 

H 

enthalpy  flow  (J  s-1) 

Hb 

upper  heating  value  flow  (J  s-1) 

i 

current  density  (mA  cm-2) 

io 

exchange  current  density  (mA  cm-2) 

I 

electric  current  (A) 

k 

thermal  transmission  coefficient  (W  K-1  m-2) 

l 

characteristic  length  (m) 

m 

mass  flow  (g  s-1) 

M 

molar  mass  (g  mol-1) 

n 

moles  of  exchanged  electrons  per  mole  of  the 
fuel  (mol  mol-1) 

Nu 

Nusselt  number  (-) 

P 

pressure  (bar) 

P 

power  (W) 

Pr 

Prandtl  number  (-) 

Q 

heat  flow  (J  s-1) 

r 

evaporation  heat  (Jkg-1) 

R 

specific  ohmic  resistance  (£2  cm2) 

Re 

Reynolds  number  (-) 

t 

temperature  (°C) 

T 

temperature  (K) 

U 

voltage  (V) 

V 

volume  flow  (lmin-1) 

w 

velocity  (ms-1) 

X 

absolute  humidity  (kg  kg-1) 

Greek  symbols 

a 

heat  transfer  coefficient  (W  m-2  K-1) 

p 

volume  expansion  coefficient  (-) 

£ 

emission  ratio  (-) 

0 

position  of  the  surface  (-) 

X 

heat  conduction  resistance  (W  m-1  K-1) 

V 

kinematic  viscosity  of  the  flowing  substance 
(m2  s-1) 

V 

process  gas  utilisation  (%) 

P 

density  (kg  m-3) 

X 

time  (s) 

<p 

relative  humidity  (-) 

Subscript 

amb 

ambient 

air 

air 

ce 

cell 

conv 

convection 

cons 

consumption 

cw 

cooling  water 

FC 

fuel  cell 

fore 

forced  convection 

free 

free  convection 

h2 

hydrogen 

in 

inlet 

loss 

losses 

max 

maximum 

out 

outlet 

rad 

radiation 

rev 

reversible 

s 

saturation 

St 

steam 

obtained  simulation  results.  Such  a  simplified  model  will  be 
presented  in  the  following  paragraphs. 


2.  Test  equipment 

The  experimental  application  for  the  investigations  can  be 
seen  in  Fig.  1.  The  fuel  cell  test  bench  is  equipped  with  a 
programmable  logic  controller  which  overtakes  the  control 
and  safety  tasks. 

For  the  definition  of  the  measurement  algorithms  and  the 
acquisition  of  the  measurement  values  a  PC  with  the  neces¬ 
sary  data  acquisition  and  the  measurement  application  was 
installed. 

The  electronic  load  was  controlled  by  the  PC  and  a  func¬ 
tion  generator  Agilent  33120  A.  The  graphical  environment 


PC  with 
data  acquisition 


Fig.  1.  Test  application  for  the  experimental  investigations. 
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VEE  Pro®  was  used  to  define  the  measurement  algorithms. 
To  prepare  the  fuel  cell  model  a  stack  was  measured  from  the 
open  circuit  operation  to  the  point  of  rating  at  variation  of  the 
process  parameters. 

The  appropriate  measurements  like  stack  voltage,  temper¬ 
atures,  pressures,  humidifier  temperatures  as  well  as  hydro¬ 
gen  and  air  utilisation  were  acquired  over  a  time  period  of 
1  min  at  current  steps  of  1  A.  The  hydrogen  utilisation  was 
changed  between  40  and  60%  in  10%  steps,  the  air  utilisation 
from  10  to  25%  in  5%  steps,  the  stack  temperature  from  40 
to  65  °C  in  5  °C  steps  and  the  air  humidifier  temperature  was 
changed  by  2  °C  steps. 


3.  Fuel  cell  model 

3.1.  Model  of  the  U-I- characteristic 

The  characteristic  of  a  PEM-fuel  cell  can  be  subdivided 
into  three  characteristic  areas. 

If  only  the  two  technically  relevant  areas  are  considered 
the  characteristic  can  be  emulated  with  the  following  empir¬ 
ical  approach  [1,2]: 

Uce  =  U0-h\ni-  Ri  (1) 

with 

Uo  =  t/rev  +  b  In  io  (2) 

where  I/rev  is  the  reversible  cell  potential  and  io  and  b  are 
Tafel  parameters  of  the  oxygen  reaction  [3].  The  linear  area 
of  the  characteristic  is  mainly  influenced  by  the  ohmic  re¬ 
sistance  R.  This  ohmic  resistance  contains  the  ohmic  losses 
within  the  membrane,  the  activation  resistance  of  the  hydro¬ 
gen  oxidation  reaction  and  the  mass  flow  resistance  on  the 
oxygen  electrode. 

The  approach  can  be  used  if  the  cells  show  the 
same  behaviour  over  the  whole  operating  range.  The  be¬ 
haviour  can  be  assumed  for  optimised  and  well-engineered 
stacks. 

Fig.  2  shows  the  dependence  of  the  stack  voltage  on  the 
current  for  the  investigated  stack  at  different  fuel  cell  tem¬ 
peratures. 

Since  the  cells  show  a  slightly  different  behaviour  and 
the  operating  parameters  could  not  be  kept  completely  con¬ 
stant  the  different  terms  of  the  regression  approach  do  not 
exactly  reflect  the  single  over- voltages.  Furthermore,  at  dif¬ 
ferent  stack  temperatures  and  higher  air  utilisation  voltage 
drops  can  occur  since  a  definite  transport  of  the  produced 
water  cannot  be  ensured  because  of  a  low  pressure  drop.  At 
higher  temperatures  and  fitted  air  volume  flows  good  or  very 
good  results  can  be  achieved  with  the  regression  approach 
used. 


Fig.  2.  U-I  characteristic  and  determined  regression  at  different  stack  tem¬ 
peratures  (air  utilisation  17%,  hydrogen  utilisation  50%). 

3.2.  Temperatures 

The  temperatures  of  the  reactants  on  the  respective  outputs 
of  the  stack  were  measured  with  temperature  sensors  which 
belonged  to  humidity  measurement  equipment  from  the  com¬ 
pany  Testo.  These  Nil 0000  sensors  were  soldered  together 
with  the  humidity  sensor  on  a  thread  block.  This  block  was 
screwed  into  a  specially  designed  mounting  which  enabled 
the  process  gases  to  directly  pass  by  the  sensors. 

Since  the  anode  outlet,  the  outlet  for  the  cooling  water 
and  the  inlet  for  the  cooling  water  are  all  situated  on  the  same 
stack  side  it  is  useful  to  have  a  representation  of  the  outlet 
temperature  of  hydrogen  dependent  on  the  outlet  tempera¬ 
ture  of  the  cooling  water  and  of  the  outlet  temperature  of  the 
air  dependent  on  the  inlet  temperature  of  the  cooling  water 
(Figs.  3  and  4).  It  could  be  noticed  that  there  was  a  linear 
dependence  of  the  cooling  water  temperatures  and  the  outlet 
temperatures  of  the  process  gases  which  could  be  emulated 
by  a  linear  regression  approach. 

During  the  measurement  of  the  characteristics  it  could  be 
observed  that  the  outlet  temperatures  of  the  gases  increased 
with  increasing  current  density  at  constant  air  and  hydrogen 
utilisation  and  constant  stack  temperature  by  approx.  1-2  K. 
This  temperature  change  can  be  attributed  to  the  imperfect 
insulation  surrounding  the  temperature  sensors  since  there 
were  no  remarkable  temperature  changes  in  dependence  on 
air  and  hydrogen  utilisation  at  nominal  power  or  high  volume 


Fig.  3.  Hydrogen  outlet  temperature  dependent  on  the  cooling  water  outlet 
temperature  and  the  hydrogen  utilisation. 
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Fig.  4.  Air  outlet  temperature  dependent  on  the  cooling  water  outlet  tem¬ 
perature  and  the  air  utilisation. 

flow  of  the  process  media.  For  this  reason  the  outlet  tempera¬ 
tures  on  the  anode  are  lower  than  on  the  outlet  of  the  cathode 
because  of  the  lower  gas  flow. 

3.3.  Pressure  losses 

The  supplied  air  has  to  fulfil  several  tasks.  Firstly,  the 
oxygen  is  made  available  for  the  chemical  reaction  with 
the  air  transport.  The  oxygen  concentration  on  the  outlet  is 
lower  than  on  the  inlet  since  the  reactant  is  consumed  via  the 
length  of  the  cathode  channel.  Lower  oxygen  concentration 
results  in  higher  activation  over- voltages  and  diffusion  losses 
[4].  The  decrease  of  concentration  can  be  influenced  by 
the  stoichiometry  number  on  a  limited  scale.  The  influence 
of  the  air  flow  on  the  oxygen  concentration  decreases  with 
increasing  stoichiometry  number.  The  intermediate  oxygen 
concentration  in  the  air  via  the  transport  path  amounts  to 
10%  at  a  stoichiometry  number  of  1.  A  concentration  of 
15%  is  reached  at  a  stoichiometry  number  of  2  and  17.5%  is 
adequately  reached  with  a  number  of  4. 

Furthermore,  the  humidity  share  in  the  cell  is  regulated 
with  air  flow  variation.  The  ability  of  air  to  transport  water 
increases  super-proportionally  with  increasing  temperature. 
Additionally,  the  cathode  air  flow  transports  a  higher  amount 
of  water  with  decreasing  air  utilisation  and  reaches  the  point 
of  saturation  via  longer  channel  length  [5]. 

Higher  stoichiometry  numbers  at  constant  channel  cross- 
section  cause  higher  flow  velocities  and  higher  pressure 
losses.  Simultaneously,  the  transport  ability  of  the  air  vol¬ 
ume  flow  at  higher  pressures  is  lower  so  that  the  problem  of  a 
dehumidified  membrane  at  the  air  inlet  at  counter-flow  of  the 
reactants  can  be  reduced  at  a  defined  air  volume  flow.  Fur¬ 
thermore,  with  increasing  pressure  difference  it  is  easier  to 
remove  the  liquid  water,  which  can  be  found  inside  the  stack. 

However,  high  pressure  differences  demand  a  high  flow 
rate  of  the  compressor  and  therewith  a  higher  energy  demand. 

The  design  of  the  flow  field  has  a  large  influence  on  the 
behaviour  of  the  fuel  cell  characteristic  and  the  electrical  ef¬ 
ficiency  of  the  facility. 

However,  the  abilities  of  the  flow  field  design  are  also 
limited.  For  example,  channels  with  narrow  cross-sections 


Fig.  5.  Pressure  losses  on  the  anode  side  dependent  on  the  hydrogen  volume 
flow. 

can  lead  to  non-uniform  supplied  zones  or  cells  because  of  the 
larger  pressure  difference  between  inlet  and  outlet  of  the  fuel 
cell  stack  [6] .  Furthermore,  the  effectiveness  of  the  reactant 
supply  and  the  water  removal  are  not  only  dependent  on  the 
geometric  parameters  like  channel  length  and  height,  but  also 
on  the  velocity  of  the  entering  gases  [7] . 

If  the  stack  is  a  commercial  product  and  the  flow  field 
is  therewith  defined,  only  the  volume  flow  dependent  on  the 
stack  temperature  and  the  stack  power  can  be  adapted.  For  this 
case  the  dependence  of  the  pressure  losses  via  the  stack  on  the 
air  and  hydrogen  flow  has  to  be  investigated  (Figs.  5  and  6). 

On  the  anode  side  there  are  normally  lower  pressure  losses 
than  on  the  cathode  side  since  the  volume  flow  of  the  hydro¬ 
gen  is  by  far  lower  than  the  air  volume  flow. 

3.4.  Heat  losses 

The  heat  losses  between  the  stack  and  the  surrounding 
environment  are  caused  by  convection  and  radiation. 

The  heat  released  via  the  stack  surface  by  convection  can 
be  calculated  with  Eq.  (3)  if  the  heat  transfer  coefficient  of 
the  stack  surface  and  the  temperature  difference  between  sur¬ 
roundings  and  surface  are  known  [8]: 

0FC,conv  —  OtAft'pC  —  Lmb)  (3) 

The  convective  heat  transfer  and  therewith  the  heat  transfer 
coefficient  are  influenced  by  a  number  of  physical  properties 
of  the  flowing  medium.  The  most  important  parameters  are 


Fig.  6.  Pressure  losses  on  the  cathode  side  dependent  on  the  air  volume  flow. 
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the  density,  the  specific  heat  capacity  at  constant  pressure, 
the  heat  conductivity,  the  volume  coefficient  of  expansion, 
the  kinematic  viscosity  and  the  thermal  diffusivity. 

The  heat  transfer  coefficient  at  free  convection  can  be  cal¬ 
culated  from  the  Nusselt  number,  the  heat  conductivity  of  the 
flowing  medium  as  well  as  the  characteristic  dimensions  of 
the  body  [9] : 


NuX 

a  =  - 

/ 


The  flow  velocity  is  a  result  of  the  lifting  power  at  free  flow 
and  the  Nusselt  number  can  be  calculated  from  the  Grashof- 
and  the  Prandtl-number: 


Nu  =  0C(Gr  Pr)" 


a 

a> 

*  i-H 

o 
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a> 
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Fig.  7.  Calculated  heat  transfer  coefficient  at  free  and  forced  convection 
dependent  on  the  stack  temperature  (ambient  temperature  25  °C). 


The  constants  C  and  n  depend  on  GrPr  and  6  considers  the 
location  of  the  surface.  The  Prandtl-  and  the  Grashof-number 
are  defined  in  the  Eqs.  (6)  and  (7): 

v 

Pr  =  -  (6) 

a 


Gr  = 


FC  Gmb)P 


According  to  Zhukauskas  the  Nusselt  number  on  horizontal 
tubes  can  be  calculated  with  Eqs.  (8)  and  (9)  for  the  convective 
heat  transfer  and  at  free  flow  [10]: 


Nu  =  0.28  Re 0  6  Pr036 


Re  = 


wd 


v 


Re  >  1000 


(8) 

(9) 


If  one  presents  the  heat  transfer  coefficient  for  free  and  forced 
convection  dependent  on  the  temperature  with  the  assump¬ 
tion  that  the  heat  transfer  shows  a  similar  behaviour  to  the 
convection  on  a  cross  passed  tube  it  can  be  shown  that  the 
heat  transfer  coefficient  for  a  forced  convection  shows  only  a 
slight  change  within  the  investigated  temperature  range  and 
the  coefficient  for  the  free  convection  only  slightly  increases 
from  a  temperature  difference  of  approx.  4-5  K. 

It  can  also  be  assumed  that  the  heat  losses  caused  by  con¬ 
vection  increase  nearly  proportionally  with  the  increasing 
temperature  difference  (Fig.  7). 

The  stack  used  in  this  study  contains  neither  a  cross-flow 
against  a  tube  nor  an  ideally  passed  body. 

Furthermore,  the  stack  is  connected  in  a  heat  conducting 
manner  via  the  mounting  and  the  media  and  current  support¬ 
ing  tubes  to  the  surroundings. 

An  experimental  investigation  of  the  heat  losses  via  the 
stack  surface  is  necessary. 

The  heat  losses  which  are  caused  by  the  heat  radiation  can 
be  calculated  with  the  following  approach: 


0FC,rad  =  Apc^FC  ,  am  bCs 


(10) 


where  £Bz,amb  is  lhe  resulting  emission  ratio  between  the  fuel 
cell  surface  and  the  surfaces  of  the  surroundings. 


The  losses  of  heat  radiation  are  nearly  linearly  dependant 
on  the  stack  temperature  (Fig.  8)  and  have  nearly  the  same 
value  as  the  convection  losses  if  a  constant  temperature  of 
the  surrounding  surfaces  and  low  fuel  cell  temperatures  can 
be  assumed. 

In  principle  there  are  two  possible  experimental  ap¬ 
proaches  for  evaluating  the  heat  losses  via  the  fuel  cell  sur¬ 
face.  In  the  first  case  the  stack  is  heated  and  kept  at  a  defined 
temperature  while  the  gas  supply  is  switched  off  and  the  re¬ 
leased  heat  flow  is  measured  via  the  stack  inlet  and  outlet 
temperatures  as  well  as  the  volume  flow  of  the  heating  water. 
This  method  was  only  minimally  suitable  since  the  temper¬ 
ature  difference  of  the  cooling  flow  on  the  inlet  and  outlet 
of  the  stack  was  very  low.  This  low  temperature  difference 
resulted  in  strongly  changing  values  for  the  calculated  heat 
flow. 

For  getting  some  reliable  results  the  heat  flow  was  mea¬ 
sured  10  times  for  15  min  at  defined  stack  temperatures  and 
the  mean  value  and  the  standard  deviation  were  determined 
afterwards. 

The  second  and  better  possibility  is  to  record  the  cooling 
down  curve  and  calculate  the  heat  losses  with  the  recorded 
temperature  curve  and  the  heat  capacity  of  the  stack. 

For  these  investigations  the  stack  was  heated  up  to  a  tem¬ 
perature  of  70  °C  and  after  the  heat  supply  the  temperature  of 


Temperature  Difference  (K) 

Fig.  8.  Calculated  heat  loss  flow  caused  by  radiation  and  convection  losses 
dependent  on  the  stack  temperature  (ambient  temperature  25  °C,  emission 
coefficient  0.76,  heat  transfer  coefficient  5.2  Wm-2  K). 
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Temperature  Difference  (K) 

Fig.  9.  Heat  loss  flow  dependent  on  the  temperature  difference  determined 
from  the  cooling  down  curve. 
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Fig.  10.  Energy  flows  of  the  stack. 


the  stack  was  measured  at  defined  time  intervals.  The  ambi¬ 
ent  temperature  was  approx.  23  °C  and  only  slightly  increased 
during  the  experiment  to  24  °C.  The  heat  capacity  of  the  fuel 
cell  could  be  defined  with  0.68  kg  kg-1  K_1  based  on  the 
known  volume  shares  of  the  single  stack  components  and  the 
respective  specific  heat  capacities  and  densities. 

Starting  with  an  assumed  linear  temperature  decrease  in 
the  time  intervals  the  heat  losses  can  be  evaluated  dependent 
on  the  change  of  the  temperature  difference  between  stack 
and  surroundings  and  related  to  the  time  interval: 


GfC,1oss 


mFCTFcApFC  lamb) 

Ar 


(11) 


As  the  calculations  for  the  heat  losses  of  the  stack  separated 
into  convection  and  heat  radiation  denote  there  is  an  almost 
linear  dependence  of  the  heat  loss  flow  on  the  temperature 
difference  (Fig.  9). 

The  comparison  of  the  heat  losses  acquired  from  the  cool¬ 
ing  curve  of  the  stack  and  calculated  from  the  heat  flow  trans¬ 
ported  by  the  cooling  water  shows  a  correlation  of  the  values 
of  both  methods  [11]. 

The  corresponding  regression  line  can  then  be  calculated 
with  Eq.  (12): 


GfC,1oss  — 


-2.571  +  1.51  |  fcw.FC,°ut  +  ^cw.FC.in  _ 


(12) 


3.5.  Energy  balance 


calculated  with  Eq.  (14): 


^H2,FC,cons 


IMu2 
2  F 


(14) 


The  enthalpy  flow  coming  from  the  outlet  of  the  anode  can 
be  calculated  via  the  mass  balance: 


#H2,FC,out  =  (fi*H2,FC,in  “  ^H2,FC,cons) 

X  [Cp,H2,FC,ouEH2,FC,out 

+  Th2,FC  ,ou  tO"0,FC  ,out 

+  Cp, St, FC, outfit, FC, out)]  (15) 

Beside  the  sensitive  heat  the  evaporated  heat  is  also  consid¬ 
ered  in  the  enthalpy  flow. 

The  demand  of  chemical  energy  can  be  calculated  via  the 
theoretical  cell  voltage  (voltage  related  to  the  upper  heating 
value)  or  the  demanded  hydrogen  mass  flow  and  the  upper 
heating  value: 

Hb  =  fi*H2,FC,cons^0  (16) 

The  provided  air  which  is  transported  through  the  humidifier 
can  be  considered  nearly  saturated. 

The  enthalpy  of  the  gas-steam  mix  results  from  the  sum 
of  the  enthalpies  of  the  mixture  ingredients: 

7^air,FC,in  =  Fair,  FC,  in  Pair,  in  x  [Cp,air,FC,in^air,FC,in 

+  Ts,FC,in(^0,FC,in  +  Cp,St,FC,in4ir,FC,in)]  (17) 


For  the  energy  balance  a  total  of  eight  energy  fluxes  have 
to  be  recorded  which  are  represented  in  Fig.  10. 

The  hydrogen  is  supplied  to  the  stack  in  nearly  dry  state 
if  the  humidifier  is  switched  off.  The  corresponding  enthalpy 
flow  can  be  calculated  with  Eq.  (13)  if  the  standard  volume 
flow,  the  temperature  and  the  pressure  on  the  inlet  of  the  anode 
are  measured: 

#H2,FC,in  =  EH2,FC,in  PH2,FC,in  Cp,H2,FC,in  fli2,FC,in  (13) 

The  mass  flow  of  consumed  hydrogen  is  directly  proportional 
to  the  current  if  the  gas  crossover  is  neglected  and  can  be 


The  parameters  of  the  air  and  steam  like  density  and  heat 
capacity  can  be  calculated  from  regression  approaches  which 
are  available  in  literature  sources. 

The  vapour  pressure  which  is  necessary  for  determining 
the  absolute  humidity  and  is  dependant  on  the  temperature 
can  be  also  calculated  from  a  regression  approach  or  by  using 
the  Antoine-equation. 

With  the  relative  humidity  and  the  vapour  pressure  the 
share  of  humidity  in  the  air  can  be  determined: 

X  =  0.622  —  (18) 


P  ~  VPs 
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The  enthalpy  flow  at  the  cathode  outlet  can  now  be  calculated 
since  the  amount  of  water  supplied  via  fresh  air  as  well  as  the 
amount  of  water  produced  by  the  chemical  reaction  is  known 
and  the  humidity  levels  at  the  outlet  of  the  electrodes  have 
been  measured.  The  water  removed  from  the  electrode  is  not 
completely  gaseous.  For  this  reason  a  term  has  to  be  provided 
in  the  equation  which  considers  the  liquid  share  of  water  in 
the  waste  air: 

^air,FC,out  =  771air,FC,out[F/?,air,FC,oudair,FC,out 

T  -^air,FC,out(^0,FC,out  T"  <%,St,FC,out^air,FC,out) 

T  (-^air,FC,max  —  -^air,FC,out)^w^air,FC,out]  (19) 

The  specific  heat  capacity  of  the  waste  air  is  now  calculated 
from  the  mixture  ratio  of  oxygen  and  nitrogen. 

This  mixture  ratio  is  also  considered  for  the  calculation  of 
the  absolute  humidity  since  the  specific  gas  constant  of  the 
waste  air  is  slightly  smaller  than  the  constant  of  the  fresh  air. 
Beside  the  heat  losses  the  final  parameter  of  the  power  has  to 
be  calculated  for  the  energy  balance: 

Pfc  =  Upchc  (20) 

The  heat  flow  of  the  cooling  water  result  from  Eq.  (21): 

Gcw,FC  =  H"  E^H2,FC,in  H-  77air,FC,in  —  ^H2,FC,out 

•  _  • 

-  #air,FC,out  ~  ^FC  “  C?FC,loss  (21) 


Fig.  12.  Calculated  and  measured  heat  flow  via  the  cooling  water  (stack 
temperature  55  °C,  hydrogen  utilisation  50%,  air  utilisation  14%,  air  inlet 
temperature  22  °C). 

3.6.  Comparison  of  model  and  realistic  behaviour 

The  applicability  of  the  simplifications  still  has  to  be 
proven  since  a  reduced  set  of  parameters  with  single  values 
was  used. 

Figs.  13  and  14  show  the  calculated  cooling  heat  flow  and 
the  regression  of  the  measured  cooling  heat  flow  for  both 
presented  cases. 

The  tolerance  of  both  heat  flows  ranges  from  approx.  2  to 
5  up  to  200%  whereby  the  absolute  failure  does  not  amount 
to  more  than  10  W. 


The  exemplary  results  of  the  measured  and  calculated  cool¬ 
ing  water  flow  for  the  carried  out  measurements  are  presented 
in  Figs.  11  and  12  for  a  stack  temperature  of  approx.  55  and 
45  °C,  a  hydrogen  utilisation  of  50%  as  well  as  an  air  utilisa¬ 
tion  of  17  and  14%,  respectively,  and  an  air  inlet  temperature 
of  22  and  26  °C,  respectively. 

The  dew  point  of  the  fresh  air  is  also  situated  at  these 
temperatures.  The  calculated  cooling  heat  flow  is  based  on 
the  aforementioned  approaches.  The  measured  cooling  heat 
flow  was  interpolated  with  a  regression  approach  of  second 
order.  The  graphs  of  the  calculated  and  measured  cooling  heat 
flow  are  nearly  congruent  and  support  the  proposed  approach 
as  well  as  the  metering  precision  of  the  single  sensors. 


Fig.  13.  Modelled  heat  flow  via  the  cooling  water  based  on  the  data  set  from 
Table  1  and  heat  flow  interpolated  from  measured  data. 


Stack  Current  (A) 

Fig.  11.  Calculated  and  measured  heat  flow  via  the  cooling  water  (stack 
temperature  45  °C,  hydrogen  utilisation  50%,  air  utilisation  17%,  air  inlet 
temperature  26  °C). 
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Fig.  14.  Relative  deviation  of  measured  and  calculated  heat  flow  via  the 
cooling  water. 
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Table  1 


Reduced  parameter  set 


Stack  temperature  (°C) 

45 

55 

Hydrogen  inlet  temperature  fH2,FC,in  (°Q 

25 

25 

Hydrogen  inlet  humidity  <^H2,FC,in 

0 

0 

Hydrogen  outlet  humidity  <PH2,FC,out 

0.82 

0.3 

Hydrogen  utilisation  vho  (%) 

50 

50 

Air  inlet  temperature  fair, fc, in  (°C) 

22 

26 

Air  inlet  humidity  <pair,FC,in 

1 

1 

Air  outlet  humidity  <pair,FC,out 

0.82 

0.78 

Air  utilisation  vair  (%) 

17 

14 

Regression  parameter  for  the  characteristic 

Cell  voltage  Uce 

958.7 

981.2 

Tafel  parameter  b  (mV) 

33.277 

34.157 

Internal  resistance  R  (£2  cm2) 

0.2780 

0.3039 

The  reduced  parameter  set  for  the  fuel  cell  investigated  in 
the  paper  is  presented  in  Table  1 . 

4.  Conclusion 

If  the  operation  of  a  PEM-fuel  cell  is  considered  under 
stationary  conditions  with  the  characteristics  of  a  fuel  cell,  the 
energy  balance  can  be  described  with  a  relatively  low  number 
of  data.  The  heat  losses  on  the  stack  surface,  the  pressure 
losses,  the  stack  outlet  temperatures  and  the  characteristic  can 
be  simulated  by  using  regression  approaches.  The  use  of  these 
regression  approaches  result  in  a  reduced  set  of  parameters. 

The  inlet  temperatures,  the  humidity  levels  of  the  pro¬ 
cess  gases,  the  stack  temperature  as  well  as  the  volume  flows 
and  the  humidity  levels  on  the  outlet  of  the  stack  form  the 


set  of  necessary  parameters.  With  the  mentioned  regression 
approaches  and  parameters  the  energy  flows  can  now  be  cal¬ 
culated.  The  results  reflect  the  energy  balance  and  the  energy 
flows  on  the  interfaces  to  the  surroundings  in  an  adequate 
manner.  This  can  be  used  in  combination  with  energy  mod¬ 
els  of  the  peripheral  components. 
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